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We present the first modeled search for gravitational waves using the complete binary black hole gravita- 
tional waveform from inspiral through the merger and ringdown for binaries with negligible component spin. 
We searched approximately 2 years of LIGO data taken between November 2005 and September 2007 for sys- 
tems with component masses of 1-99 Mrr^ and total masses of 25-100 Mg. We did not detect any plausible 
gravitational-wave signals but we do place upper limits on the merger rate of binary black holes as a function of 
the component masses in this range. We constrain the rate of mergers for l9Ma, < mi , 1112 < 28Mp;binary black 
hole systems with negligible spin to be no more than 2.0 Mpc^-' Myr^'at 90% confidence. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.60.Jd, 97.60.Lf, 97.80.-d 



I. INTRODUCTION 

This paper presents a search for gravitational waves 
from binary black hole (BBH) coalescences with total 
mass 25 Mp, < M < 100 Mp, and component masses IMq < 
m\,m2 < 99 Mq. The search used Laser Interferometer 
Gravitational-wave Observatory (LIGO) [IJ data taken during 
the fifth science run (S5) from November, 2005 to Septem- 
ber 2007 when both LIGO sites were operating. The first 
LIGO site in Hanford, Washington hosts two interferomet- 
ric gravitational-wave detectors, a 4km detector, HI, and a 
2km detector, H2. The second site in Livingston, Louisiana 
hosts a single 4 km detector, LI. The Virgo detector |2| in 
Cascina, Italy commenced its first science run (VSRl) on May 
18, 2007 and since then LIGO and Virgo have operated their 
instruments as a network. However, this search did not use 
Virgo data, because it was not as sensitive to these high mass 
systems during VSRl. Additionally the GEO600 detector in 
Germany was functioning during S5. However, GEO600 data 
was not analyzed for similar reasons. The search results for 
compact binaries with total mass M < 35Mq in LIGO S5 and 
Virgo VSRl data have been reported in (3^-5]. To date no 
gravitational waves have been directly observed. 

The gravitational-wave driven evolution of BBHs is con- 
ventionally split into three stages - inspiral, merger and 
ringdown (IMR). The gravitational- wave signal during the 
adiabatic inspiral phase can be described by post-Newtonian 
(PN) expansion. This technique is very accurate for 
comparable-mass systems at large separations, but breaks 
down near the innermost stable circular orbit (ISCO) lll21l . 
Modeling of the merger requires the numerical solution of 
the full Einstein equations in a highly dynamical strong-field 
regime. After the merger, the rapidly damped quasi-normal 
ringdown of the black hole (BH) toward a stationary Kerr 
black hole is described by black-hole perturbation theory. 
This is the first analysis that incorporates a template family 
of waveforms modeling all three stages of BBH coalescence. 
This search covers systems for which the effects of BH spins 
can be neglected for detection. For BHs in this mass range 
the merger occurs in the LIGO detectors' most sensitive fre- 
quency band. 



A. Motivation to searcli for higlier mass systems 

Black holes observed in X-ray binaries range up to ^ 2OM0 
lI^HSl, and predictions from population-synthesis models have 
typically suggested the masses of components of BH-BH bi- 



naries that merge within lOGyr will mostly lie in the range 
5^0^'«i,'«2^ IOM0 Homtl. However, a number of chan- 
nels have been suggested through which significantly more 
massive black-hole binaries could form. 

Observations of ICIO X-1, a binary with a massive black- 
hole (> 24 Mp ) accreting from a Wolf-Rayet companion star, 
and a similar recently-observed binary NGC 300 X-1 |[T2l , 
suggest that more massive BH-BH binaries can form through 
isolated binary evolution, with component masses ~ 2OM0 
ifTSll . Meanwhile, several simulations over the past few years 
have indicated that dynamical formation could significantly 
contribute to coalescence rates involving BH-BH binaries in 
dense stellar environments, such as globular and nuclear star 
clusters lfT4UT7l . The most massive black holes are likely to 
sink to the centers of clusters through mass segregation and 
substitute into binaries during three-body encounters, thus fa- 
voring relatively massive components in dynamically formed 
BH-BH binaries. Moreover, the BH merger products in such 
dense clusters can be reused if they are not ejected from the 
cluster due to recoil kicks, leading to higher-mass mergers in 
subsequent generations; component masses for BH-BH merg- 
ers in globular clusters can therefore range to ^ 30Mq ifTsll 
and beyond. Additionally, although stellar winds in high- 
metallicity environments may prevent the formation of mas- 
sive black holes, mass loss through stellar winds would be 
much less significant in low-metallicity environments, allow- 
ing more massive black holes to form lITSllTQi . 

Binaries including an intermediate-mass black hole 
(IMBH) having a mass 5OM0 < m < 5OOM0 could represent 
another exciting source for LIGO and Virgo detectors. Obser- 
vational evidence for IMBHs is still under debate (see reviews 
Il20l l2ri for additional details), although a recently discovered 
ultra-luminous X-ray source |22| represents a possible IMBH 
detection. If IMBHs do exist in the centers of some globular 
clusters, they could contribute to coalescence rates in one of 
three ways: (i) through inspirals of stellar-mass compact ob- 
jects into IMBHs |l23l|24l; (ii) through mergers of two IMBHs 
that formed in the same stellar cluster 1251 : and (iii) through 
mergers of IMBHs from two different globular clusters when 
their host clusters merge |26|. It may also be possible to de- 
tect mergers of binary IMBHs arising from a direct collapse 
of early population III binary stars |27|. Although the rates of 
events involving IMBHs are highly uncertain, they may reach 
tens of detections per year in the Advanced LIGO/Virgo era 
(see 1 28 1 for a review of expected detection rates for all bi- 
nary types relevant to the advanced-detector era). 

Advances in numerical relativity (NR) in the last 5 years 
have enabled this search for BBH coalescence to significantly 
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improve on previous methods and results. We describe those 
advances and the impact on gravitational-wave astronomy in 
the next section. 



B. Complete inspiral, merger and ringdown waveforms from 
numerical relativity 

Constructing nonperturbative numerical solutions for the 
merger of two black holes has proven to be remarkably dif- 
ficult. It has taken more than four decades since Hahn and 
Lindquist first attempted the numerical investigation of col- 
liding black holes |291 to compute the gravitational-wave sig- 
nal from the last orbits, merger and ringdown of a black- 
hole binary system. These simulations are now possible us- 
ing many different methods. Only months after Pretorius' ini- 
tial breakthrough in 2005 ll30l . his success was repeated us- 
ing a different approach lISTl l32l . and since then several NR 
groups were able to produce increasingly accurate BBH sim- 
ulations exploring increasing portions of the parameter space 
(see, e.g. Il33l436l for recent overviews on the field, and Sec. 
2 of f37l for a compact summary). The success of NR sim- 
ulations has lead to a range of new physical insights includ- 
ing the calculation of recoil velocities produced by asymmet- 
ric emission of gravitational radiation during the merger pro- 
cess 1 38 - 56] and the prediction of the parameters of the rem- 
nant Kerr BH for a wide class of initial configurations P0ll57l - 
l62l 1621 469 II. Most importantly for the gravitational wave (GW) 
community, these simulations were able to compute accurate 
GW waveforms from the late inspiral and merger of BBHs 
in many configurations, and these predictions were success- 
fully matched to PN and effective one-body (EOB) predic- 
tions |70-88|. Thus, by combining analytical and numerical 
calculations, it is now possible to construct accurate waveform 
templates coherently modeling the IMR of the coalescence of 
binary black holes (BBHs), as described in section [II] 



C. Summary of past searches 

The LIGO and Virgo Scientific collaborations previously 
searched for systems that are a subset of the BBH parame- 
ter space explored by this analysis using different techniques. 
Previous searches used templates that modeled only the in- 
spiral or ringdown phases. None employed templates that in- 
clude IMR waveforms. To date no gravitational waves from 
BBH coalescence have been observed. 

The first search explored systems with component masses 
in the range SM© < mi, m2 < 20Mq in ^ 386 hours of LIGO's 
second science (S2) run [89] with a 90% sensitivity to sys- 
tems up to 1 Mpc. The second search covered sources up 
to 40 Mq total mass for LIGO's third science run (S3, 788 
hours) and up to 80 Mq in total mass for the fourth science 
run (S4, 576 hours) [90| with ^ 10 times the range sensitiv- 
ity of S2, and placed a 90% confidence upper limit on the 
merger rate of ^ 0.3Lj"q^ yr^^ for systems with a total mass 
of ~ 4OM0. (Here Lio is lO'" times the blue Solar luminos- 
ity and is used as a proxy for the expected number of sources 



in a galaxy. For searches that extend beyond ^ 20 Mpc there 
are approximately 0.0198Lio Mpc^^|91 1. The presently dis- 
cussed search extends beyond 20 Mpc and has sufficient sen- 
sitivity to use units of Mpc^^ Myr^'. Additionally the Milky 
Way is 1.7Lio.) The S2, S3 and S4 science runs used 
phenomenological waveforms proposed by ll92l that extended 
the inspiral to higher frequency but did not include the com- 
plete IMR signal nor the effects of spin. A search that used 
templates including the spin effects was conducted over S3 
data targeting asymmetric systems with component masses 
IMq < nil < 3M0 and 12M0 < 1112 < 20Mq |93|. 

BBH mergers with sufficiently high mass will have most 
of their in-band gravitational-wave amplitude in the ringdown 
phase. A search over S4 data probed the ringdown phase of 
BBH coalescence and placed 90% confidence upper limits on 
the merger rate for systems with total masses 85 Mq < M < 
39OM0 of 1.6 X 10--^Lio' yr-i (32Mpc-3 Myr"') lH- The 
ringdown waveforms are a function of the final state of the 
black hole and do not depend on the details of the merger 
For this reason a spectrum of initial states (arbitrary spins 
and component masses) can be probed via a ringdown-only 
search. 

Finally, the LIGO Scientific Collaboration (LSC) and Virgo 
searched S5 data and Virgo's first science run data (VSRl, 
which overlapped with the last ~ 6 months of S5) for BBHs 
with a total mass up to 35 Mq The 90% confidence 

upper limit on merger rate for black hole binaries with total 
mass ^ 30 was ^ 3 x IO^'^Ljq' yr"\ which is ^ 6 Mpc"-' 
Myr-i 0. 



D. Summary of the present results 



No plausible gravitational-wave signals were detected in 
this search. The loudest events are discussed in section |IV] 
Despite not detecting BBH signals directly we are able to 
infer an upper limit on the merger rate of such systems in 
the nearby Universe. We do not impose a particular popu- 
lation model within our mass range and instead present our 
merger rate limits as a function of component mass ranges. To 
90% confidence we constrain the rate of mergers for 19M0 < 
mi, 1112 < 28M0binary black hole systems to be no more than 
2.0Mpc^^ Myr^'. We highlight numbers from this range be- 
cause it may include some of the heavier BH binaries that may 
arise in population synthesis models (e.g., LJ9 I), but was not 
covered by the S5 low-mass search lHJQ. Additional mass 
pair rate limits are given in section [V] 

The paper organization is as follows. In section |II] we de- 
scribe two families of waveforms used in this search (effective 
one-body model and phenomenological IMR model). Section 
III summarizes the key points of our data analysis pipeline. 
This part includes information about the template bank, data 
quality, background estimation and candidate ranking statis- 
tics. Section IV contains results of the search, in particular we 



present the loudest events. In section|V]we discuss detection 
efficiency, which we estimate by injecting simulated signals 
into the detector data. Additionally in this section, we present 
an upper limit on the coalescence rate for this search. Finally, 
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section VI presents the conclusions and plans for future im- 
provements. 



n. WAVEFORMS USED IN THIS SEARCH 



Modeled waveforms are invaluable tools for extracting 
weak signals from noisy data in gravitational-wave searches 
for compact binaries |95|. The models are used to efficiently 
filter the data for signals and to assess the sensitivity of the 
instruments and data analysis procedure via simulations. This 
section motivates the need for new waveform models in this 
search and describes the models chosen. 

For binaries with a total mass in the range targeted by this 
search, 25 — IQQMq, the ISCO is reached in the sensitive fre- 
quency band of the LIGO interferometers. Thus, standard 
inspiral-only PN waveforms, which are typically terminated at 
the ISCO frequency, do not capture all of the observable sig- 
nal. Furthermore, this abrupt, in-band end of the search tem- 
plates can be problematic for the signal consistency checks. 
On the other hand, IMR templates model all of the observable 
signal and naturally decay away during the ringdown phase 
rather than abruptly ending. For these reasons, it is highly de- 
sirable to use IMR template waveforms to search for binary 
coalescences in this mass range. In figure [T] we plot an ex- 
ample IMR waveform in the time and frequency domains and 
note the extra signal relative to inspiral-only waveforms. 

Fortunately, the recent breakthroughs in numerical relativ- 



ity (see Sec. IB I have revealed the nature of the merger and 
ringdown phases of BBH coalescences. While it is infeasi- 
ble to use the NR simulations directly as search templates, in- 
sights gained from the simulations have informed the develop- 
ment of analytic IMR waveform models. Currently, two main 
paradigms exist in the construction of IMR waveforms. In the 
BOB approach, an effective-one-body description of the two- 
body problem is tuned with NR simulations and then matched 
to the quasi-normal modes of the BH ringdowns to produce 
analytical IMR waveforms in the time domain. In the phe- 
nomenological IMR model, the NR waveforms are matched to 
PN waveforms to produce "hybrid" PN-NR waveforms which 
are then parametrized to produce analytical IMR templates 
in the frequency domain. The BOB waveforms are used as 
search templates and also as injected waveforms to test our de- 
tection efficiency. The phenomenological waveforms are used 
for injections and provide a check that our search pipeline can 
detect waveforms which are slightly different than our search 
templates. The next two subsections describe each of these 
families of analytic IMR waveforms. 

Since BOB waveforms are generated in the time domain, 
the presence of an abrupt starting point at a given low fre- 
quency can result spurious high frequency power. To mitigate 
these effects a tapering window was applied to the beginning 
of the generated BOB waveforms [96J . 
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FIG. 1: Example of the EOBNRvl IMR waveforms used in this 
search for a (25 + 25) Mr., binary optimally located and oriented at 
100 Mpc in the time domain (top panel) and the frequency domain 
(bottom panel). The solid vertical lines mark the location of the 
Schwarzschild ISCO, which is the termination point for inspiral-only 
waveforms. 



A. Effective one-body model 

The BOB approach, originally introduced in ||97ll98l , pro- 
vides a PN-resummed Hamiltonian which can be used to 
evolve a binary system through its inspiral and the final 
"plunge" of the compact objects before they merge. This 
trajectory can be used to generate a waveform (f) 
which can be matched onto a waveform /I'^rger-RD^f^ describ- 
ing the merger and ringdown of the resulting black hole, made 
up of a superposition of the black hole's quasi-normal modes. 
The two pieces are combined at a suitably chosen matching 
time fmatch to produce an inspiral-plunge-merger-ringdown 
BOB waveform I 



meiaer-RD 



0(f- 



'match 



(11. 1) 
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where 0() is the Heaviside step function. 

The inspiral-plunge EOB waveform at leading-order ampU- 
tude in a PN expansion is determined from the trajectory r(f ), 
0(f) as [98] 



/jinsp-plunge^j-j 



4GM77 (GMdi^VI^ 



(11.2) 



where is the luminosity distance. We now summarize the 
fundamentals of the EOB calculation of the trajectory; more 
details can be found in [79-86, 97-100|. As usual, my andm2 
are the black hole masses, M = mi+ m2 is the total mass of 
the binary, jj. — mi 1712/ M is the reduced mass, and rj — ji /M 
is the symmetric mass ratio. 

For a binary with negligible spin effects, the motion is con- 
fined to a plane and can be described in the center of mass by 
polar coordinates (r, 0). The conservative dynamics is then 
captured by a Hamiltonian H^'^^{r,pr,p,j,). The trajectory is 
evolved according to Hamilton's equations 



dr _ ^//EOB 
dt d Pr 
d^ _ ^//EOB 

dt ^ 

dpr 

dt 
dp^ 
dt 



{r,Pr,P(l>), 

{npr,p,i,), 



dr 

^(bir.Pr.Pt,) 



{r,Pr,P^) 



(II. 3a) 

(II.3b) 

(11.3c) 
(II.3d) 



The inspiral of the binary comes about due to the addition of 
non-conservative dynamics in the last of Hamilton's equations 
via the tangential radiation-reaction force arising from the 
basic PN expression of the energy flux. Here we use a Ke- 
plerian Pade resummation IIOIJ of the energy flux as given 
by Eq. (15) of ISTj . More recent models have used more so- 
phisticated fluxes, such as the p-resummation [102] and non- 
Keplerian flux models which describe non-quasi-circular ef- 
fects 1 82, 83 , 85 , 86 1. 

The form of the EOB (resummed) Hamiltonian is ||97ll 



(11.4) 



where H^^ is the effective Hamiltonian 

A{r) pI , 



1 



pI 



+ 2(4-3tj)tj 
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(11.5) 



and where the radial potential functions A (r) and D{r) appear 
in the effective metric ll97l 



dslff = -A(r) dt^ + dr^ + (dO^ + sin^ d(j)^ 



Air) 



The Taylor-approximants to the coefficients A (r) and D(r) can 
be written as 



k+i 

i=0 
k 



GMY 
rc^ I 



OlM = I*(l7) 



rc^ J 



(II.7a) 
(II.7b) 



The functions A(r), D{r), Ai.(r) and /^^.(r) all depend on 
the symmetric mass ratio rj through the T7-dependent coef- 



2GM 



and 



ficients a; (77) and £/,(?]). [When T] — > 0, A(r) — > 1 — 
D{r) —i' 1 and the metric (11.61 reduces to the Schwarzschild 



metric] These coefficients are currently known through 3PN 
order (i.e., up to A; = 3) and can be found in [81 [. Dur- 
ing the last stages of inspiral and plunge, the EOB dynam- 
ics can be adjusted closer to the numerical simulations by in- 
cluding in the radial potential A(r) a pseudo 4PN coefficient 
05(77) = 05 T], with 05 a constant. Here, we follow [81 [ and fix 
05 — 60. We refer to this model, the first NR-adjusted EOB 
model implemented for a search of GW data, as EOBNRvl. 
Since [^Tl was published, more accurate numerical simu- 
lations became available and more sophisticated EOB mod- 
els have been calibrated. This includes a different value of 
a$ j|82l484l [86^ and also the introduction of a pseudo 5PN co- 
efficient 06(77) = ae 77 [85 [, with ag a constant. We refer to 
the second NR-adjusted EOB model implemented for a search 
of GW data, as E0BNRv2. This most recent EOB template 
family has been developed in II103II : it includes the latest im- 
provements [82 -84 86] to the EOB model and also other re- 
finements which are necessary to match highly-accurate NR 
waveforms for a broad range of mass ratios. 

In order to assure the presence of a horizon in the effective 
metric (II. 6 1, a zero needs to be factored out from A (r). This 
is obtained by applying a Pade resummation t 99| . The Pade 
resummations of A(r) and D{r) at pseudo 4PN order are de- 
noted A4(r) and 1)4 (r) II122I . and the explicit form used in this 
paper can be read from lISTl . 

The merger-ringdown waveform in the EOB approach is 
built as a superposition of quasi-normal modes, [79, 981 I104II 



n=0 



^match ) 



(11.8) 



where n is the overtone number of the Kerr quasi-normal 
mode, is the number of overtones included in our model, 
and A„ are complex amplitudes to be determined by a match- 
ing procedure described below. We define (J„ = o,, — /a„, 
where the oscillation frequencies a)„ > and the inverse 
decay-times a„ > 0, are numbers associated with each quasi- 
normal mode. The complex frequencies are known functions, 
uniquely determined by the final black-hole mass and spin. 
They can be found in 01O5L The final black-hole masses and 
spins are obtained from the fitting to numerical results worked 
out in [81 [. 



The complex amplitudes A„ in Eq. (II. 8 1 are determined by 



(11.6) 



matching the EOB merger-ringdown waveform with the EOB 
inspiral-plunge waveform close to the EOB light ring. In par- 
ticular, here we use the matching point which is provided an- 
alytically by Eq. (37) of 181 1 . In order to do this, we need 
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independent complex equations that are obtained at the match- 
ing time by imposing continuity of the waveform and its time 
derivatives. 



(11.9) 



with A: = 0, 1,2, • • • jA' — 1. In this paper we use N=3. The 
above matching approach is referred to as point matching. Al- 
though it gives better smoothness around the matching time, 
it is not very stable numerically when is large and higher- 
order numerical derivatives are employed. More sophisticated 
matching procedures have been proposed in the literature to 
overcome the stability issue 1821 l83l l85l [86l . and will be 
adopted in the future. 



B. Phenomenological IMR model 



Reference II 10611 presented a different way of constructing 
non-spinning IMR waveforms by combining PN calculations 
with numerical simulations. They first constructed a family 
of hybrid waveforms by matching PN waveforms with NR 
waveforms in certain overlapping time intervals where both 
the approaches are expected to be valid |87|. Restricted 3.5PN 
waveforms in the TaylorTl approximation were matched to 
NR waveforms produced by the BAM NR code |107|. These 
hybrid waveforms were used to construct a family of analyt- 
ical waveforms in the Fourier domain, of the form h{f) = 
'^eff(/)s''^^^'''^'i where the effective amplitude and phase are 
expressed as: 



[ (// /merg) 



-7/6 



/</r 



merg 



[ Wif (/, /ring, C7) /ring < / < /cut 

^eff(/)^ ^ t(^kri'+ykn+zk){nMff-'y' 



+ 27tfto + (pQ. 



(n.IOb) 



In the above expressions, C is a numerical constant whose 
value depends on the location and orientation of the binary as 
well as the physical parameters, -Sf (/,/ring, c) is a Lorentzian 
function that has a width a, and that is centered around 
the frequency /ing. The normalization constant w is cho- 
sen so as to make Aeff(/) continuous across the "transition" 
frequency /-ing. The parameter /merg is the frequency at 
which the power-law changes from /^^/^ to f^'^^^. The phe- 
nomenological pai-ameters Hj = {/merg, /ring, O'j/cut} are given 
in terms of the physical parameters of the binary as: nMpLj = 
ajT]^ +bjri +Cj. The coefficients {aj,bj,Cj\j = 0. . .3} and 
{^k,yk,Zk\k = 0,2,3,4,6,7} are tabulated in Table I of II108I . 
We refer to the waveform family defined by these coefficients 
as IMRPhenomA, and these are the waveforms used for injec- 
tions in the present search. These waveforms are generated in 
the frequency domain and are then converted to the time do- 
main for injections by means of the inverse Fourier transform. 



The choice of the time interval for matching PN and NR 
waveforms is somewhat ad hoc. Currently, the matching in- 
terval is chosen so as to maximize the fit of PN and NR wave- 
forms. Moreover, the PN waveforms employed in the match- 
ing are computed in the restricted PN approximation, and the 
amplitude of the NR waveforms is scaled to match with PN 
waveforms. This causes the amplitude of the waveforms to 
have a systematic bias of ^ 10%. Later improvements in this 
model have already addressed some of these issues ifSSl : we 
refer to this improved waveform family as IMRPhenomB . 



C. Systematic errors in waveform models 

Although the two waveform families have been tested via 
comparisons to numerical waveforms, there are a number of 
possible sources of systematic uncertainty in the two wave- 
form families. For example, as discussed above, there are 
subtleties in choosing the matching interval between PN inspi- 
ral waveforms and numerical simulations when constructing 
the hybrid waveforms used to calibrate the Phenomenologi- 
cal waveforms. Similarly, the BOB procedure to attach the 
merger-ringdown waveform to the inspiral-plunge waveform 
can be quite delicate and become unstable if not done prop- 
erly. Most notably, both waveform families have been tested 
against numerical simulations only in the nearly-equal-mass 
regime, up to ratios of 3 : 1 and 4:1. It is not clear whether 
these waveforms are faithful to the actual signals in the case 
of highly unequal masses. 

The waveform models initially used for this search, EOB- 
NRvl and IMRPhenomA, had both been revised by the time 
the search neared completion. The revisions, which in- 
cluded improved cahbration and more accurate matching to 
NR waveforms, as well as improved modeling of the post- 
Newtonian inspiral phase, resulted in waveforms that were 
significantly more faithful. For example, the revised version 
of the phenomenological waveform family, IMRPhenomB, 
has systematic biases in signal-to-noise ratio (SNR) of < 10% 
relative to NR waveforms in the mass range of interest and for 
mass ratios below 4:1; for the revised version of the FOB 
waveform family, E0BNRv2, the systematic biases in SNR 
relative to NR waveforms are < 3%. These systematic bi- 
ases do not account for any errors in the NR waveforms them- 
selves. 

The largest effect of the revision of both models was to 
systematically reduce the gravitational-wave amplitude dur- 
ing merger We found that within our errors it was sufficient 
to adjust the distance of the simulated signals to take into ac- 
count the lower intrinsic gravitational-wave amplitude in the 
corrected models. The upper limits quoted in section [V| are 
thus based on a search carried out with EOBNRvl templates, 
but with the distances of EOBNRvl and IMRPhenomA injec- 
tions adjusted to match the SNR of the revised E0BNRv2 and 
IMRPhenomB waveform models. 

We can get a sense of the systematic uncertainty in the 
waveform amplitudes by comparing the SNRs between the 
two waveform families. We find that the SNR of the most 
recent versions of the two families, E0BNRv2 and IMRPhe- 
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noniB, agrees to better than ^ 10% for mass ratios less than 
6 ; 1 in the mass range of interest, but diverges by nearly 50% 
for mass ratios of 10 : 1. The latter value is chosen as the limit 
on the mass ratio for phenomenological IMR injections. 

III. THE DATA ANALYSIS PROCEDURE 

The data analysis procedure involves a multi-stage pipeline 
that automates the extraction of signals from the data, the anal- 
ysis of coincident events between detectors and the estimation 
of background. The pipeline used for this search was similar 
to that of previous S5 searches fJ'-'Sl except for the choice of 
template waveforms and some minor points described in sub- 
sequent sections 

Here we summarize the data analysis procedure. First, data 
for the three different detectors, HI, H2, and LI are divided 
into 2048 s blocks in order to estimate the time dependent 
power spectral density (PSD) of the detector noise. The PSD 
is required to choose the search templates and to filter the data 
itself. Next, the data are processed in a two-stage procedure. 
The first stage filters the data with the templates and identifies 
potential events in each detector. Then the pipeline checks for 
coincidence between detectors. We allow double and triple 
coincident combinations between detectors. After finding co- 
incident events the data are re-filtered using only the templates 
that participated in the coincident events. The data needs to 
be filtered with fewer templates at the second stage after de- 
manding coincidence with other detectors. The second filter- 
ing stage employs the veto L109J, which drastically reduces 
the background of this search, but is too computationally ex- 
pensive to be performed during the first filtering stage with the 
full template bank. Once coincident events are identified they 
are clustered in a ten second window to produce a maximum 
of one coincident event every ten seconds. We apply the same 
procedure to time-shifted data streams and compare the time- 
shifted results to the zero-lag results to assess the significance 
of our events. This procedure is repeated with simulated sig- 
nals in order to assess the sensitivity of the pipeline. 

In the remaining sections, we elaborate on this procedure 
emphasizing differences with previously published searches. 

A. Generation of coincident event candidates 

In this section we describe the process of obtaining candi- 
date events. First we discuss how to choose templates to filter 
the data. Next we describe the filter process itself and how 
to identify events that are significant in a single detector. We 
then describe how we check for coincident events between 
detectors. Finally we describe how data quality impacts our 
assessment of candidates. 



1. Selection of search template parameters 

The observed gravitational waveform depends on the com- 
ponent masses of the binary. A bank of template waveforms 



called a template bank is chosen to adequately cover the pa- 
rameter space of possible waveforms. The template bank 
used for the search consisted of templates covering total mass 
between 25-100 Mq, and component masses between 1-99 
Mq. The bank was tiled using a hexagonal placement algo- 
rithm 111 lOL such that the intended minimum SNR was 97% of 
its maximal value [1 1 1 1. The template spacing was determined 
using the metric calculated for the stationary phase approx- 
imation 111 121 extended to the effective ringdown frequency. 
This metric, terminated at ISCO frequency, was used in pre- 
vious searches for signals from low-mass systems. Although 
the metric is not formally correct for the FOB templates used 
in this search it has been found that the bank provides the de- 
sired minimal match for most of the parameter space and at 
worst a 95% match for the high-mass region of the bank. The 
average number of templates required to cover this parameter 
space was ^ 1600 per detector 

2. Filtering 

After properly selecting the templates to cover the mass pa- 
rameter space, the data are filtered. The signal to noise ratio 
for a given template waveform h{t) is a convolution of the 
template with the data weighted by the noise power spectral 
density, defined as 

z{t)^4r^^e'''-f'df, (in.la) 
^0 S„[f) 

^0 S„{f) 

where the tilde and * denote a Fourier transform and a com- 
plex conjugate, s refers to the data and Sn is the noise power 
spectral density. In order to cover the entire parameter space 
all of the data are filtered with every template. z{t) is a com- 
plex time series where the real part represents a template phase 
of and the imaginary part represents a phase of n/2. The 
real-valued SNR p is given by p{t) = \z{t)\/a. 

We trigger on the local maxima of each filter's time series 
when the SNR is above 5.5 and record those times and tem- 
plate parameters. A list of triggers is then passed to the next 
stage of the pipeline, which checks for coincidence between 
detectors. 



3. Coincidence test 

We require events to be coincident in at least two detectors. 
For events to be considered coincident, the time of coales- 
cence and the masses of the svstem lll23l given by the trig- 
gers in each instrument must agree to within a certain toler- 
ance |1 13J . Although we allowed for double coincident com- 
binations, we discarded H2L1 events that lacked an HI trigger 
if HI was operating. Since HI was more sensitive than H2 it 
should have produced a trigger for a real event. 

As with the searches for low mass compact binary coales- 
cences in S5, we used a coincidence test based on the template 
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bank metric. This test accounts for correlations between the 
different parameters and attains a lower false alarm rate for a 
given detection efficiency than simple parameter cuts. As was 
noted earlier, the metric used in this search was suboptimal. 
To take into account this limitation, the coincidence require- 
ments were looser than those of previous S5 searches ISJQ. 



4. Data quality vetoes 



Detectors Analyzed time (yr) 

H1H2L1 0.6184 

HlLl 0.0968 

H2L1 0.0609 

TABLE I: The analyzed time surviving the pipeline after category 3 
vetoes were applied. H1H2 times were not analyzed due to the in- 
ability to properly estimate the background for co-located detectors. 



Not all of the data taken during S5 was used for this anal- 
ysis. The detectors frequently lost lock or were taken out of 
lock for commissioning work. Only times with stable lock 
stretches deemed as analyzable were marked as science time. 
Segments of science time containing more that 2048 s of data 
were analyzed in each of the three detectors HI, H2, LI. 

Occasionally data quality (DQ) during science time suf- 
fered from transient excess noise. Significant work was done 
to characterize these times prior to examining the search can- 
didates so as to not bias our detection and upper limit state- 
ments |114|. Events at times suffering from poor data qual- 
ity are removed from the analysis. The procedure of vetoing 
events reduces the live time and also the false alarm rate of the 
search |114|. The following describes the basic procedure for 
vetoing candidates based on DQ. 

The detectors are sensitive to a variety of noise transients 
(glitches) of non-astrophysical origin, such as instrumental 
glitches and environmental disturbances. The status of the de- 
tectors is monitored by a number of auxiliary data channels 
that record the internal degrees of freedom of the interferome- 
ters and the output from environmental sensors. When the sta- 
tus of a detector is suboptimal, the time is flagged. Because 
the templates used in this search have an impulse response 
lasting ^ 10 seconds, a short glitch can produce triggers last- 
ing several seconds after the glitch occurs. DQ flag intervals 
often require search-specific time padding to improve the ef- 
fectiveness of the flag. The length of this padding is deter- 
mined by looking at the distribution of triggers in the flagged 
interval. The effectiveness of a DQ flag is evaluated by the 
following metrics: Efficiency: the percentage of single de- 
tector triggers flagged. Because these triggers are analyzed 
before coincidence, they are dominated by transient noises lo- 
cal to the detector; Dead-time: the percentage of flagged time; 
Used-percentage: the fraction of flags that contain at least one 
background trigger An effective flag has a high efficiency, a 
high used-percentage and a low dead-time. Flags found to be 
effective by these metrics are used as vetoes. 

DQ flags are classified into four veto categories accord- 
ing to their metrics' performance. Category 1 contains times 
when the data was not analyzed at all as described at the be- 
ginning of this section. Category 2 includes vetoes with a 
high efficiency-to-dead time ratio and a high used percentage. 
The origin of these glitches is well-understood and time in- 
tervals are well-tuned. Category 3 vetoes times with noise 
sources whose coupling with the gravitational-wave channel 
is less understood, such as those due to environmental noise. 
Category 3 vetoes are less correlated with transients and are 
characterized by higher dead time and lower used percentage 



than category 2 vetoes. Some flags, for example the overflows 
of digital channels monitoring the alignment of the interfer- 
ometer arm lengths and mirrors, belong to both category 2 
and 3 with different window lengths. Category 4 contains ve- 
toes with low efficiency and high dead time. These flags usu- 
ally identify minor environmental disturbances and problems 
recorded in the electronic logbooks. 

The DQ vetoes are used in the following way. Category 

2 vetoes are used unconditionally in the search. We exam- 
ine events after Category 2 for detection candidates. However 
we apply Category 3 vetoes before creating the list of candi- 
dates used to constrain the BBH merger rate. The category 

3 veto list is chosen in advance in order to not bias our rate 
limit results. Category 4 vetoes are used only to follow up in- 
teresting candidates, they do not have any impact on the rate 
limits quoted in this paper All Category 2 or greater vetoes 
are applied after the second coincidence stage before cluster- 
ing to produce the event list. Vetoed time is accounted for to 
ensure that the analyzed time calculations are correct. Table |l] 
gives the analyzed time available after Category 3 vetoes are 
applied. 



B. Ranking and evaluation of candidate events 

1. Signal consistency check 

Accounting for data quality as described in the previous 
section is not sufficient to remove all triggers caused by en- 
vironmental or instrumental sources. For that reason we em- 
ploy a two stage pipeline that performs an additional signal 
consistency check. In the second filtering stage we explicitly 
check the match of the signal to the template waveform by per- 
forming a chi-squared test II109I . In this test, the template is 
divided into p frequency bins (for this search, we use 10 bins) 
such that each bin contains the same expected contribution to 
the total SNR, if the signal matches the template exactly. The 
SNR of the trigger in each bin is compared to the expected 
SNR, and the differences are added in quadrature to obtain 
the value of x^- We decompose the template waveforms into 
p pieces of identical power (J^/p 

1=1 

Uiif) = h{f) &{f - /,■ low) ©(/;,high - /) . (in.2b) 
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Using (III. la I we compute a filter time series for each of the 
orthogonal pieces. 



(III.3) 



The X statistic is then computed as 



P 



(111.4) 



Since z(f )is a complex number, corresponding to both phases 
of the filter, the statistic has 2{p — 1) degrees of freedom. 

Previous searches in this mass range did not use IMR wave- 
forms. Since the models were not accurate they did not use a 
test I.90J . The statistic already provides significant sep- 
aration from noise for a large fraction of simulated signals 
in this search. Future search efforts in this mass range might 
employ new signal-based vetoes and multivariate classifiers to 
achieve a better separation of signal from background lll24l . 

Once the statistic is evaluated we have almost all of the 
inform ation ne cessary to begin ranking events. We describe in 
how the X^ statistic is folded together with the 



section 



111B3 



SNR to produce a ranking statistic known as effective SNR. 
First, however, in section jlll B 2| we describe how we estimate 
our background, which is also required for ranking the coin- 
cident events. 



2. Background estimation 

We assume that instrumental noise triggers are not cor- 
related between detectors. We estimate the background of 
this search by examining accidental coincidences from time 
shifted data. This section describes how we estimate the back- 
ground. The next section describes how the background esti- 
mate is used in ranking events. 

In order to estimate the background of coincident events we 
repeat the coincidence analysis with 100 time shifts between 
the two LIGO sites in multiples of five seconds. We call the 
events found by this procedure time slide events. We expect 
that there will be no correlated noise between the sites. There- 
fore, the time-shifted analysis provides 100 background trials 
to which we then compare the un-shifted data. Unfortunately 
the assumption of uncorrected noise was not adequate for the 
collocated Hanford detectors, HI and H2. All events found in 
HI and H2 but not LI were discarded due to correlated noise 
corrupting the background estimate. 

We find that the estimated background of the search is a 
function of time, the parameters of the signals searched for 
and which detectors observed the event. The total mass of the 
recovered signal is the best single parameter that tracked the 
signal parameter dependence of the time slides. We elaborate 
in sectionii 
events. 



Ill B 3 how this was used in the ranking of candidate 



Observation 


Detectors 


Detectors 


Mass 


Epoch 


Functioning, F Participating, P 


Range, M 


12 ~two-month 


HlLl 


HlLl 


[25,5O)M0, 


epochs 


H2L1 


H2L1 


[5O,85)M0, 




H1H2L1 


H1H2L1 


and 




HlLl 


[85,1OO]M0 



TABLE II: Breakdown of the analysis time and coincident trigger 
sets. The LIGO S5 run is divided into 12 epochs, each roughly two 
months in duration; within each epoch, the time is divided accord- 
ing to which detectors were operating and not vetoed. Since there 
are three allowed combinations of functioning detectors, there are 
12 X 3 = 36 different analyzed time periods. Different combinations 
of coincident events are allowed depending on which detectors are 
functioning and participated in the coincident event. There are a to- 
tal of four possible functioning/participating detector combinations 
which contribute to the analysis. Within each observation epoch 
and functioning/participating detector combination, the events are di- 
vided into three mass bins according to the average total mass of the 
templates involved in the coincident event. This means there are a 
total of 12 X 4 X 3 = 144 different types of coincident events. Each 
type of event has a separate background distribution used to calculate 
its false-alarm probability. 



3. Ranking events 

The ranking of candidate events is a multi-stage process. 
The end ranking statistic is a false alarm rate (FAR) for each 
event that indicates how often events like it (or louder than it) 
occur in time slides. This section describes how we compute 
the FAR and rank our events. 

First, single-detector triggers are assigned an effective SNR 
Peff which is a function of p and x^- The functional form is 
chosen to match the false alarm rate contours of the single- 
detector background in the SNR - x^ plane. The effective 
SNR is defined as 



Peff = 



[(i+pV50)(zVzif)] 



1/4 ' 



(III.5) 



where 50 is an empirically determined parameter and xlof = 
2{p—\) is 18 for this search. The single detector effective 
SNRs Peff.j- are combined in quadrature to give a coincident 
effective SNR 



Peff.c 



EPeff.i 



(III.6) 



We compute the FAR by comparing the un-shifted events to 
the time slide events. Due to the non-Gaussian properties of 
the detector noise, the FAR depends on the template. It also 
depends on how many detectors were operating and partici- 
pated in the event. We compute the FAR as a discrete func- 
tion of four parameters, the total mass M, the detectors that 
participate in the coincidence P, the detectors that were func- 
tioning but not vetoed at the time of the coincidence F, and the 
combined effective SNR rank of the event R. We will denote 
a time-slide event that estimates our background as B. Each 
parameter is an index for the event B. The first and second 



13 



indices, F and P, describe the instraments that were function- 
ing during the event and the detectors that participated in the 
event. Only the following combinations were considered: 1) 
triggers found in HI and LI when only the HI and LI de- 
tectors were operating 2) triggers found in HI and LI when 
all three detectors HI, H2 and LI were operating 3) triggers 
found in H2 and LI when only H2 and LI were operating and 
4) triggers found in all three detectors when all three detectors 
were operating. Note that as mentioned previously we were 
not able to estimate a reliable background for triggers found 
only in HI and H2. Therefore those events were discarded. 
We also discarded events found in H2 and LI when all three 
detectors were on since the more sensitive HI should observe 
a real signal. To summarize, the following shorthand notation 
for the 4 combinations of participating P and functioning F 
detectors will be used: P,F e {H1L1,H1L1; H1L1,H1H2L1; 
H2L1,H2L1; H1H2L1,H1H2L1}. The third index M denotes 
a range for the total mass estimated for the event and is in 
the set {[25,50), [50,85), [85, 1OO)}M0. The fourth index 7? is 
the rank of the event given by its effective SNR, Peff,c- The 
R index is determined by assigning the event having a given 
P, F and M with the lowest combined effective SNR defined 
in ( |ni.6| l the value and the the next lowest 1, etc., until all 
events are ranked. We calculate the false alarm rate FAR for 
a given event as the number of all time slide events, B, with 
a rank (/?+) larger than that event's rank divided by the time 
analyzed Tf in the time shifted analyses, which is a only a 
function of the instruments that were on and not vetoed. 



FAR 



PFMR 



^ BpFMR+ Tf 



(III.7) 



R+>R 



This now allows us to map a zero-lag (unshifted) event to a 
FAR by assigning it the same four parameters . 

In addition to the indices describing how the false alarm rate 
was computed, there is one remaining implicit parameter that 
refers to the time of the events. We separated the two calen- 
dar years of data into 12 two-month periods. Each was treated 
separately for the calculation of ( |I11.7| i in order to crudely cap- 
ture the variation of the noise properties over the course of 
S5. It is worth making explicit the number of combinations 
over which false alarm rates were computed. Each of the 12 
two-month periods had 4 possible combinations of detectors 
that were functioning and that produced triggers as mentioned 
above. Additionally each had 3 total mass bins. The result is 
12x4x3 = 144 separate calculations of ( |111.7| i. This is de- 
scribed additionally in table |ll] and is relevant for interpreting 
the significance of events in section [TV| 

Next we assess the FAR of the events independently of 
the mass range M and the participating detectors P in order 
to compute a global ranking that only takes into account the 
detectors that were functioning and no other parameters. To 
do this we use the inverse FAR^' as an intermediate ranking 
statistic to replace combined effective SNR as the rank in the 
index R. We denote these newly ranked time slide events as 
B'. Then the combined FAR is 



The combined FAR is only a function of the detectors that 
were functioning F during the event and the inverse FAR rank 
computed at the previous step R. From the combined FAR 
we can also compute the False Alarm Probability (FAP). As- 
suming Poisson statistics, we define the FAP as the chance 
of getting one or more events louder than the event in ques- 
tion purely from background. This is defined as FAPf^ — 
1 — exp (— FAR/trT/t). Tp is nominally a particular detector 
combination live time for a two-month analysis period, but 
can be replaced with the entire observation time in order to 
obtain the FAP for an event given the result of all 12 two- 
month periods. 



IV. LOUDEST COINCIDENT EVENTS 

As previously mentioned, we divided the ^ 2 calendar 
years of data into 12 two-month blocks and this resulted in 
144 separate computations of ( |111.7| i. Combining the FAR us- 
ing ( |111.8| ) resulted in 36 separate periods consisting of distinct 
times when a given set of detectors were functioning and pro- 
viding data. These categories are independent since they arise 
from distinct times. Since a lower FAR implies a more signif- 
icant event, we use FAR^' to rank the events. 



Table 111 gives the top 10 loudest events of the search ranked 
by FAR^^ Three of the 10 candidates were louder than any 
events in the 100 time-shifted coincident sets used to estimate 
the background. The table provides the bound on the FAR 
based on the total observed background time during their two- 
month period when the same detectors were functioning. We 
note that it is not surprising to have events louder than the 
background given the limitations of the background estima- 
tion. We used only 100 time shifts and the number of tri- 
als examined for the computation of the FAR was 144. We 
therefore expected to observe ^1.4 events more significant 
than our estimated background and we observed three. In or- 
der to estimate the significance of these three events we em- 
ployed two additional techniques. As the primary method, we 
first interpolated and extrapolated the FAR from our 100 time- 
shift background estimate. To obtain an alternative estimate, 
we extended our time-shift study to 1000 shifts for the two- 
month periods in which those events occurred. We decided 
before un-blinding the analysis to use the extrapolated FAR 
values in the upper limit computation when necessary. We 
also examined many properties of these events in a qualitative 
follow-up procedure. The result of our analysis is that all three 
events have FAPs of > 10^^, assuming the full 0.8 yr obser- 
vation time, and all are consistent with rare instrumental noise 
fluctuations; none are plausible candidate gravitational-wave 
detections. This section provides some additional detail about 
these events. 



A. H1H2L1 event at GPS time 848905672.3369 
(November 30, 2006 07:07:38.3369 GMT ) 



FAR 



FR - 



R+>R P M 



PFMR+ Tp 
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The loudest event of this search at GPS time 
848905672.3369 was found in all three detectors and 
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Rank 


FAR (yr 


') Peff 


GPS Time 


Phi 




mm 


mim 


Ph2 


Xm 






Pli 


xti 






1 


<0.20 


12.8 


848905672.3369 


172.0 


4057.9 94.0 


6.0 


24.4 


167.4 


49.7 


17.3 


8.3 


46.0 


95.2 


4.8 


2 


<0.25 


11.6 


825664840.1523 


5.6 


21.5 


51.7 


1.1 


6.2 


1.6 


50.5 


1.1 


5.5 


39.1 


36.2 


2.4 


3 


< 1.40 


10.3 


842749918.8057 


- 






- 


5.5 


7.8 


67.1 


2.5 


12.2 


20.4 


83.2 


16.8 


4 


2.7 


12.0 


830222610.4062 


5.5 


34.2 


98.0 


2.0 




* 


* 


* 


28.8 


43.5 


91.7 


8.3 


5 


5.4 


9.8 


849056023.4121 


5.7 


11.3 


29.9 


1.3 


* 


* 


* 


* 


5.6 


2.5 


23.6 


1.8 


6 


9.0 


9.8 


827865922.1265 


9.2 


67.6 


37.3 


1.2 


- 


- 


- 


- 


6.8 


4.2 


31.0 


1.5 


7 


12 


9.5 


836048263.0366 


6.2 


15.0 


52.9 


1.4 


64 


14.6 


46.7 


1.6 


5.9 


21.3 


53.0 


1.3 


8 


12 


10.7 


854487078.6543 


6.1 


29.6 


96.7 


3.3 




* 


* 


* 


18.1 


29.8 


97.0 


3.0 


9 


13 


10.8 


835998008.6890 


23.2 


52.3 


94.8 


5.2 


* 


* 


* 


* 


5.8 


21.2 


78.1 


1.2 


10 


15 


9.8 


857817894.5767 


8.8 


29.7 


90.7 


1.4 


9.9 


40.8 


94.8 


5.2 


5.9 


28.1 


90.4 


1.4 



TABLE III: The loudest events of the search. The coincident events are ranked by their combined false alarm rate FAR. A " - " represents 
that the detector was not functioning during the time of the event in question. A " * " represents that the detector was functioning but did not 
produce a trigger above the single detector SNR threshold of 5.5. Notice that the top three events were found above their local background 
estimates. For that reason only limits on their combined FARs are given here. See the text for details. 



was more significant than any of the time-shifted events in its 
background estimate. We put a bound on its FAR from the 
original 100 time-shift background estimate of 1 per 5 years. 
We also estimated the FAR by interpolating and extrapolating 
the original 100 time-shift background estimate using a fit to 
the trigger distribution. The extrapolated FAR was 1/1 .6 yr 
Note that it was larger than the bound due to the fitting 
procedure in the tail of the trigger distribution. Also note that 
we decided in advance to use the extrapolated FARs for the 
rate limit calculations in the next section. Therefore the FAR 
used for this event was 1/1.6 yr We computed 1000 additional 
time shifts in this two-month period to better estimate its false 
alarm probability. From the additional time-shift background 
estimate we computed that this event had a false alarm rate of 
1/50 yr. 

Given that we searched nearly 1 year of data, this event is 
consistent with fluctuations. The conservative probability of 
getting this event in background (by choosing the lowest of 
the FAR estimates) is ^ 0.02. Our assessment of this candi- 
date is that it is a loud glitch in HI with a moderate response 
in H2 coincident with low amplitude noise in LI. The ratio 
of distance estimates associated with the signals in HI and 
H2 II125I is not consistent with a signal. We measured a ratio 
of ~ 10 and it should be ^ 1. The HI does not lie within 
the expected signal distribution. We therefore conclude that 
this is not a gravitational- wave detection candidate. 



B. H1H2L1 event at GPS time 825664840.1523 
(March 06, 2006 07:20:26.1523 GMT) 

The second loudest event of this search at GPS time 
825664840.1523 was more significant than any of the 100 
time slides performed during the two month period in triple 
coincident H1,H2,L1 time. The event was found in all three 
detectors HI, H2, and LI with SNR only slightly above 
threshold 5.60, 6.17 and 5.55 respectively. The masses were 
consistent between the detectors. In HI and LI this event had 
a that was consistent with both time slide events and sig- 



nals. However, it had an unusually low %-■ value (0.1 per de- 
gree of freedom) in H2. Kx^ value of less than 0. 1 per degree 
of freedom is rare for both signals and noise. No background 
events out of ^ 300,000 had such a low value nor did any 
of the ^ 10^ simulated signals. The ranking of this event was 
artificially elevated by the unusually low x^ value. If this 
event had a higher of 1 per degree of freedom it would 
not stand above background. We conclude that this event is 
not a gravitational-wave detection candidate. 

The unusually low x^ value put this event in a region of pa- 
rameter space where the FAR extrapolation is not valid. This 
event happened to occur in a segment of time that we reserved 
in advance as a test data set, called a playground, that was not 
used in the rate limit calculation shown in the next section. 
See[V]for more details. We place a bound on its FAR of 1 per 
4 years from the original 100 time-shift background estimate. 
We found that this candidate is not stable to small changes in 
our analysis pipeline. We were thus not able to measure its 
FAR independently using more time slides. 



C. H2L1 event at GPS time 842749918.8057 
(September 20, 2006 01:11:44.8057 GMT) 

The third loudest event of this search at GPS time 
842749918.8057 was found in H2 and LI. It was louder than 
any of the time slide events in its two-month period. We put 
a bound on its FAR of 1.4/yr from our original analysis. An 
independent check using additional time slides yielded a FAR 
of 2.9/ yr We also interpolated and extrapolated the original 
100 time-shift run to obtain a FAR of 1.9/yr. 

The FAR of 1.9/yr was used in the upper limit calculation 
described in the next section. The LI SNR and x^ is consis- 
tent with the background in that instrument. The H2 trigger 
is just above the SNR threshold of 5.5. This event is not rare. 
With a FAR of 1 .9/yr we expected to observe an event sim- 
ilar to this in our total observation time even though it was 
above background in its local two-month H2L1 observation 
time. We conclude that this event is not a gravitational-wave 
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detection candidate. 



V. MERGER RATE LIMITS 

Before examining events for detection candidates we 
agreed upon the procedure described in this section for estab- 
lishing an upper limit on the merger rate of black hole binaries 
if no detections were found. 

In order to constrain the merger rate we had to assess the 
sensitivity of the search. To test the detection sensitivity of 
our search pipeline, we injected ~ 10* signals into the detec- 
tor strain data and processed it with the same pipeline used for 
the search. Events associated with the injected signals hav- 
ing FARs less than the loudest event of the search are con- 
sidered to be found by the pipeline. We inject both EOB and 
phenomenological waveforms into the data. The injection pa- 
rameters were as follows. For both waveform families, the 
injected signals had distances between 1 Mpc and 750 Mpc 
distributed uniformly in the logarithm of distance. Both fami- 
lies had a uniform distribution of sky location and orientation. 
For both families the total mass of the binary systems var- 
ied between 25 — IQOMq. The component mass distributions, 
however, did differ between the EOB and phenomenological 
waveforms. The component mass distribution for EOB sig- 
nals was generated by first producing a uniform distribution 
in the component masses between 1 — 99 M© and then clip- 
ping the distribution to have no systems outside of the total 
mass range 25 — 100 Mr. The mass distribution for the phe- 
nomenological waveforms was produced by first generating a 
distribution that was uniform in mass ratio (mi : OT2 ; mi > 1112) 
between 1 : 1 and 10:1 and then clipping the result to have no 
systems outside of the total mass range 25 — lOOMp,. 

As previously stated we divided the ^ 2 years of data into 
12 two-month periods and examined each of the three func- 
tioning detector combinations H1H2L1, HI LI, H2L1 sepa- 
rately for a total of 36 periods. We reserved 10% of the de- 
tector time as an un-blinded playground: we do not use play- 
ground data in computing the upper limit on the merger rate. 
The second loudest event described in section llVB] happened 



to occur in the playground time. Using (111.8 1 we ranked each 
candidate event in the 36 periods. We used the loudest event 
in the foreground after category 3 vetoes in each period to es- 
tablish a combined FAR threshold for determining what injec- 
tions were found. For the events louder than background we 
used the extrapolated FAR as agreed on prior to un-blinding 
the analysis. 

The efficiency e of recovering simulated signals in the de- 
tection pipeline is a function of the loudest event FAR, FAR*, 
the radial distance to the source r and the masses m\,m2. Note 
that in practice the mass dependence is captured by binning 
the mass plane into the boxes illustrated in figure [2] The bar 
denotes that the efficiency is averaged over sky position and 
orientation. We define the efficiency as 



e(FAR*,r,mi,m2) 



A^/-(FAR*,r,mi,m2) 



(V.l) 



M(FAR*,r,mi,m2) 
where Nf is the number of found injections, A^, is the total 



number of injections and FAR* is the FAR of the loudest event 
in a given analysis period. We then compute the volume of 
the sky surveyed in each of the 36 independent observation 
periods (denoted by the index /) by 

y,(mi,m2,FAR*) = J 47: r^Bi {FAR* ,r, mi, m2)dr, (V.l) 

which has units of Mpc^. We estimate the variance 

(j2(mi,m2,FAR*) = (y,(mi,m2)2) - (y,-(mi,m2))2 (V.3) 

by bootstrapping the input injection distribution to account for 
Monte-Carlo errors as well as varying the injection distances 
according to the conservative quadrature sum of the calibra- 
tion uncertainty among the three detectors, 20% 1 1 15| . An 
additional systematic error is associated with uncertainty in 
the target waveforms. These limits are presented with our 
best understanding of the currently available waveforms. If 
we take the fractional difference in the SNRs of Phenomeno- 
logical IMR and EOB waveforms, ^ 10% (see section 11 C I, as 



an indication of the uncertainty in the range due to imperfectly 
known waveforms, we conclude that the rates as reported in 
figure |2]have an additional systematic uncertainty of ^ 30%. 
This uncertainty is not included in the rate estimates nor are 
any other systematic errors, for example the accuracy of the 
waveform phasing. Some errors are discussed in t81iil06J . 

In order to establish a merger rate R{mi,m2) in units of 
mergers Mpc^^ yr^' we adopt formula (24) in fl 16]. It is 
important to note that some simplification of these formulas 
occurs when choosing the FAR as the ranking statistic 1 1 17 1. 
Adapting the loudest event formalism described in I1II6J to 
our notation, if we constructed a posterior on R using only the 
results of a single analysis period, the marginalized hkelihood 
function would be 



p{ki,a„A,\R)oc 



1 



+ 



RQ.iAi{l + l/ki 



{\+R£li/kif'+' {\+R£li/kif.+^ 

(V.4) 



where 



n,=\^(mi,m2,FAR*)7^-, 

^ ^ r y,(mi,m2,FAR*) 1^ 
' [(J,(mi,m2,FAR*)J ' 
^ ^ a'ln[y,(mi,m2,FAR*)] 1 
d¥AR* 



(V.5) 
(V.6) 

(V.7) 



Ti is the analyzed time for index / (assumed to have no errors), 
y is taken from ( |V.2| i, and the proportionality constant in ( |V.4[ ) 
can depend on kj and A,, but not R. 

In order to obtain the combined posterior probability distri- 
bution for the rate, given the sensitivities and loudest events 
of the 36 different analysis periods, labeled by the index we 
multiply the likelihood functions and assume an initial uni- 
form prior on the rate. This results in a posterior probability 
of the form 

p{R\mi,m2)=piR\{ki},{ni},{A,}) 

p{{ki},{ai},{Ai}\R) = llpiki,at,MR) 

(V.8) 
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FIG. 2: The 90% confidence upper limit on the merger rate as a 
function of mass in units ofMe (symmetric over mi and 1112). This 
image represents the rate limit in units of Mpc^^ Myr^' . These lim- 
its can be converted to traditional units of L^q Myr^ ' by dividing by 
0.0198Lio Mpc"^ gl). Only bins with mass ratios < 4 : 1 have up- 
per limits computed due to uncertainty in the waveform models for 
more asymmetric systems. 



FIG. 3: Average range defined in \V.9\ for the search in Mpc as a 
function of mass (symmetric over mi and 1112), assuming target wave- 
forms that match the FOB and Phenomenological models. 



the ranges found by examining the loudest event in each of the 
36 periods. 



We integrate the normalized form of ( |V.8[ ) to 90% to estab- 
lish the 90% confidence upper limit on the merger rate (still 
a function of component mass), Rgo%. The result is given in 
figure |2] The upper limit in the lowest mass bin considered 
in this search is an order of magnitude higher than the most 
optimistic binary black hole merger rates predicted by current 
population-synthesis studies (see, e.g., |[T3] [T9l |28l ). At the 
upper end of the analyzed mass range, there are no reliable 
estimates for merger rates for intermediate mass black holes, 
whose very existence remains to be confirmed; however, see 
ll24H26ll28j for some intriguing possibilities. 

As discussed above, due to the uncertainties in the wave- 
form models for asymmetric systems, we do not present up- 
per limits for mass ratios < 4 : 1 . However, we do provide an 
average range for systems with smaller mass ratios based on 
the EOB and Phenomenological waveform models, in figure 
[3] The average range is defined as 

(^(mi,m2,FAr)) = -^^7;.^,('«i,'«2,FAR*) , (V.9a) 



^,(mi,m2,FAR* 



An 



y;(mi,m2,FAR* 



1/3 



, (V.9b) 



where V,(OTi,m2,FAR*) is defined in ( |V.2| l, is the radius of 
the sphere having volume of and the average range (^) is 
the time-weighted average of ranges computed from each of 



VI. CONCLUSIONS 

We presented the result of a search for BBH coalescence 
during LIGO's fifth science run spanning approximately two 
years of data taken from fall 2005 to fall 2007. We targeted bi- 
naries with total mass M — mi +1112 in the range 25 M© <M < 
100 Mq and component masses of IMq < mi,m2 < 99M0 
with negligible spin. In order to effectually detect such sys- 
tems with LIGO it was necessary to use template waveforms 
that encompass the inspiral, merger and ringdown phases of 
compact binary coalescence. We employed two waveform 
families in this search to filter and assess the sensitivity. Both 
had been tuned to numerical relativity simulations. 

We did not detect any plausible gravitational-wave can- 
didates. However we estimated our search sensitivity and 
were able to constrain the merger rate of the targeted sources 
in the nearby Universe. We established to 90% confidence 
that the merger rate of black holes with component masses 
19M0 < mi,W2 < 28M0is less than 2.0Mpc"^ Myr"'. We 
note that this is still about an order of magnitude higher than 
optimistic estimates for such systems |28| (see also |13, 19|). 

There are a number of limitations in the current approach, 
which will be addressed in future searches. The main lim- 
itation is that the template waveforms neglect the effects of 
spin. Although the statistical distribution of the spins of black 
holes in binaries is not well known [IISJ . there are exam- 
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pies of black holes in X-ray binaries which have been ob- 
served to have a large spin 1.1 19J . For a binary with spin- 
ning components, the expected observed gravitational-wave 
signal will differ from the non-spinning case; the observed 
duration can be different and there may be modulation of the 
gravitational-wave amplitude and phase. Neglecting such ef- 
fects in the search templates will affect the detection efficiency 
for binaries with spinning components. However, at the time 
this search was conducted there were no analytical inspiral- 
merger-ringdown waveforms for systems with generic spins. 
This should be borne in mind when interpreting the results of 
the search. 

Another limitation of the search is that, due to the shorter 
duration and bandwidth of the signals in comparison to 
searches for lower mass systems, it is harder to distinguish 
between genuine signals and background events, since the sig- 
nals themselves are more "glitch-like". New approaches to the 
ranking of candidate events are being developed to improve 
the sensitivity of searches for these systems. 
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